etopo1_full = xr.open_dataset("data/ETOPO1_Bed_g_gmt4.nc") etopo1_fulletopo1_baltic = etopo1_full.sel( x=slice(9, 30), y=slice(53, 67), ) etopo1_balticetopo1_baltic.to_netcdf("data/ETOPO1_Bed_g_gmt4_baltic.nc")